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What is claimed is: 

1. A method for separating a plurality of mixed signals into a plurality of 

a 

component signals comprising the steps of: 

(a) producing a plurality of discrete Fourier transform (DFT) values corresponding 
to frequency components of an input segment of said mixed signals; 
5 (b) updating cross correlation estimation matrices using said DFT values; 

(c) computing, using a cost-function minimization process, an update value for a 
plurality of filter coefficients for a finite impulse response (FIR) filter using said cross 
correlation estimation values; 

(d) filtering said mixed signals using said FIR filter having said updated filter 
10 coefficients to separate said mixed signals into one or more component signals; and 

(e) iteratively repeating steps (a), (b), (c) and (d) for successive input segments of 
said mixed signal. 

2. The method of claim 1 wherein step (c) further comprises the substeps of: 
(cl) transforming said update filter coefficients from the frequency domain into 

the time domain; 

(c2) zeroing any filter coefficients having a value other than zero for any time that 
5 is greater than a predefined time Q, where Q is less than a value of an input-segment 
length, T; thereby producing a set of constrained time domain filter coefficients; and 

(c3) transforming said adjusted time domain filter coefficients from the time 
domain back into the frequency domain. 
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3. The method of claim 1 wherein said cost-function minimization process is a 
gradient descent process. 

4. The method of claim 1 wherein said computation of said filter coefficient 
update values further includes an adaptation of update step sizes based on a normalization 
factor. 

5. The method of claim 4 wherein said update values, including said update step- 
size adaptation, are computed according to 



where [i is a fixed learning constant, h is a weighting factor for the step-size adaptation, 
E is a filter-parameter cost function operating on a square difference respecting a diagonal 
covariance of said component signals, and is a gradient step for E. 

6. The method of claim 1 wherein a running average of said cross correlation 
estimation values are produced according to 




R x (f 9 io) = (1 - Y)R x (t,cai) + Y*(t,co)x H (f,a)) 



where x(t,co) is the mixed signal in the frequency domain. 
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7. An apparatus for separating a plurality of mixed signals into a plurality of 



component signals comprising: 



means for producing a plurality of discrete Fourier transform (DFT) values 



corresponding to frequency components of an input segment of said mixed signals; 



5 



means for updating cross correlation estimation matrices using said DFT values; 



a cost-function minimization processor for computing an update value for a 
plurality of filter coefficients for a finite impulse response (FIR) filter using said cross 
correlation estimation values; and 

an FIR filter having said updated filter coefficients to separate said mixed signals 
10 into one or more component signals. 

8. The apparatus of claim 7 wherein said cost-function minimization processor 
further comprises: 

a first transformer for transforming said update filter coefficients from the 
frequency domain into the time domain; 
5 means for zeroing any filter coefficients having a value other than zero for any 

time that is greater than a predefined time Q, where Q is less than a value of an input- 
segment length, T; thereby producing a set of constrained time domain filter coefficients; 
and 



a second transformer for transforming said adjusted time domain filter coefficients 



10 from the time domain back into the frequency domain. 
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9. The apparatus of claim 8 wherein said first transformer uses an inverse Fourier 
transform and said second transformer uses a Fourier transform. 

10. The apparatus of claim 7 wherein said cost-function minimization processor 
carries out a gradient descent process. 

11. The apparatus of claim 7 wherein said computation of said filter coefficient 
update values further includes an adaptation of update step sizes based on a normalization 
factor. 

12. The apparatus of claim 11 wherein said update values, including said update 
step-size adaptation, are computed according to 



where |i is a fixed learning constant, h is a weighting factor for the step-size adaptation, 
E is a filter-parameter cost function operating on a square difference respecting a diagonal 
covariance of said component signals, and V W E is a gradient step for E. 

13. The apparatus of claim 7 wherein a running average of said cross correlation 
estimation values are produced according to 




R x (t,co) = (1 - yyk x (t f e)) + yK{t,co)x H (t 9 &) 



where x(t,co) is the mixed signal in the frequency domain. 
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14. The apparatus of claim 7 further comprising a voice recognition system for 
processing at least one of said component signals. 

15. A computer readable storage medium containing a program that, when 
executed upon a general purpose computer system, causes said general purpose computer 
system to become a specific purpose computer system that performs a method for 
separating a plurality of mixed signals into a plurality of component signals comprising 

5 the steps of: 

(a) producing a plurality of discrete Fourier transform (DVT) values corresponding 
to frequency components of an input segment of said mixed signals; 



10 plurality of filter coefficients for a finite impulse response (FIR) filter using said cross 
correlation estimation values; 

(d) filtering said mixed signals using said FIR filter having said updated filter 
coefficients to separate said mixed signals into one or more component signals; and 

(e) iteratively repeating steps (a), (b), (c) and (d) for successive input segments of 
15 said mixed signal. 



(b) updating cross correlation estimation matrices using said DFT values; 



(c) computing, using a cost-function minimization process, an update value for a 



16. The computer readable medium of claim 15 wherein step (c) further 
comprises the substeps of: 
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(cl) transforming said update filter coefficients from the frequency domain 



into the time domain; 



5 



(c2) zeroing any filter coefficients having a value other than zero for any 



time that is greater than a predefined time Q, where Q is less than a value of an input- 
segment length, T; thereby producing a set of constrained time domain filter coefficients; 
and 

(c3) transforming said adjusted time domain filter coefficients from the 
10 time domain back into the frequency domain. 

17. The computer readable medium of claim 15 wherein said cost-function 
minimization process is a gradient descent process. ' 

18. The computer readable medium of claim 15 wherein said computation of said 
filter coefficient update values further includes an adaptation of update step sizes based 
on a normalization factor. 

19. The computer readable medium of claim 18 wherein said update values, 
including said update step-size adaptation, are computed according to 



where [i is a fixed learning constant, h is a weighting factor for the step-size adaptation, 
5 E is a filter-parameter cost function operating on a square difference respecting a diagonal 
covariance of said component signals, and V W E is a gradient step for E. 
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20. The computer readable medium of claim 15 a running average of said cross 
correlation estimation values are produced according to 

R x (t,a>) = (1 - y)R x (f 9 o)) + yK{t,a))\ H (t,co) 
where x(t,co) is the mixed signal in the frequency domain. 



